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CHAPTER I 


INTRODUCTION 

1 . 1 Previous Related Studies 

The ionosphere is the region of the upper atmosphere 
comprised of neutral and ionized gases and free electrons. 
The presence of free ions and electrons is caused mainly by 
photoionization during direct sunlight conditions, and is 
maintained during nighttime conditions by the finite recom- 
bination rates of the ionized particles. Although there are 
not any definite boundaries given to this region, it is 
generally considered to start at an altitude of 50 km and 
continue upward to the interface between the atmosphere and 
the solar wind. 

Active investigation of this region began half a 
century ago when it was observed (Appleton and Barnett 1925) 
that radio waves from ground based senders were reflected 
from an altitude of 100 to 200 km. Using this reflection 
technique, researchers were able to generate fairly detailed 
descriptions of the ionosphere. They also found that the 
ionized medium was not uniform or static, but at times con- 
tained random fluctuations in electron density which caused 
amplitude fluctuations or scintillations of the returned 
echoes. Further investigations of this effect produced some 
information on the size and motion of these inhomogeneous 
structures . 

This method contained two fundamental difficulties: 


- 2 - 

(1) it was limited to analysis of the region below the 
height of maximum ionization since no reflections could be 
obtained from greater heights, (2) the mathematical descrip- 
tion of this process near the reflection point was complex 
and obscured the interpretation of the experimental data. 

A new experimental technique was introduced by Hewish (1951) 
which eliminated the first problem and reduced the complexity 
of the second. He introduced the idea of using a radio star 
with its higher frequencies as the source in a transmission 
experiment. Further improvement has occurred in the last 
ten years with the introduction of satellite radio sources. 
These have the advantage of providing very coherent and 
stable signals with the ability to scan large sections of 
the ionosphere very quickly during the transit of the satel- 
lite across the field of view of a receiving station. 

To date most transmission experiments have dealt 
with the amplitude variations of the radio signal. Ampli- 
tude studies by Spencer (1955) showed the irregularities to 
be elongated along the magnetic field lines, and studies by 
Jones (1960) produced a range of apparent sizes of these 
structures. Further studies (Jespersen et al . , 1964; Dyson, 
1969; Kelleher et al., 1970; Rufenach et al . , 1972) have 
produced information on the statistical nature of their 
location, size, and shape. 

The experimental use of phase fluctuations rather 
than amplitude scintillation, has not received much 
attention except for a study by Porcello and Hughes (1968) 
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due probably to its much greater experimental complexity. 

As will be shown later, however, the use of phase infor- 
mation has the advantage of allowing the observation of 
larger structure sizes than by amplitude studies for a given 
source frequency and allowing a more direct interpretation 
of the results in terms of the ionospheric characteristics. 

1 • 2 Statement of the Problem 

The object of this study is to investigate the power 
spectrum of phase modulation imposed upon satellite radio 
signals by the inhomogeneous F-region of the ionosphere 
(100-500 km). Records containing the signals will be pro- 
cessed to yield temporal variations resulting predominately 
from the sweep across a receiving station as the satellite 
transits, of the two-dimensional spatial phase pattern 
produced on the ground. This will then be linked to the 
power spectrum of structure sizes comprising the diffracting 
medium. It is hoped that these results may provide useful 
information on the nature of the diffracting medium and on 
possible production mechanisms for its inhomogeneous nature. 


CHAPTER II 


THEORY OF DIFFRACTION DUE TO AN INHOMOGENEOUS IONOSPHERE 

2 . 1 Diffraction Theory 

It is necessary to develop a physical model and the 
associated mathematical description of the transmission 
experiment being conducted. The most general model is that 
of a moving point source at a finite distance from the 
diffracting medium of finite thickness. The motion of the 
source will cause the geometry of the problem to change by 
changing the angle of incidence upon the medium and the 
distances between source and diffracting medium and observer. 
The simplified model to be used is of a uniform, unit ampli- 
tude plane wave normally incident upon a thin diffracting 
medium. Later the restriction to use of a plane wave 
(source at infinite distance) will be lifted. Use will be 
made of diffraction theory (outlined by Ratcliffe 1956) to 
describe the physical model used. 

The diffracting medium modifies both the amplitude 
A and phase cf> of the plane wave, so that as it leaves the 
medium it possesses a complex amplitude f (x,y) across the 
plane parallel to the medium (complex quantities are indi- 
cated by a bar) given by: 

f(x,y) = A(x,y) exp [i<f> (x ,y) ] l\ 1 

The wave field immediately below the plane of the medium, 
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and hence at all points in this half -space, can be repre- 
sented by an angular spectrum of plane waves A C 6 ^ ® 2 ^ 
incident upon the diffracting medium at the angles 0^*02 
(6 is the complement of the angle between the propagation 
vector k and the x axis, similarly for © 2 and the y axis). 
What this involves is taking the vector sum of the angular 
components shifted in phase by an appropriate amount, with 
x =y=o taken as the reference point for zero phase shift. 
Letting C. = cos 0. and S. = sin 0. for j = 1,2, A(0,,0~) 
becomes 

A (0 1 ' © 2 ^ = ^ ^1 '^ 2 ^ 

The complex amplitude in terms of its angular spectrum is 
(lower case dimensions and coordinates x,y,z are normalized 
in terms of wave length A, upper case letters represent 
unnormalized distances) 

f(x,y) = J f F(S 1 ,S 2 ) exp[2Tti(S 1 x + S 2 y)]ds 1 ds 2 

2.2 

As the wave front leaves the medium the individual elements 
of F(S 1 ,S 2 ) are shifted in phase by an amount 2uC^z . As 
used here, the sines of © 1 and © 2 along with C 3 form the 
direction cosines of k(e l ,0 2 )« c 3 is related to 0 1 and © 2 
by 
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C3 ~ t" cos ^2 9 1 + 2 9 2^ cos *-2 6 1 ' 2 + 0 2^ 


or 


1 

c 3 = [cos(e 1 + e 2 ) cos (e ^ - e 2 )] 2 


For a given distance from the screen, equation 2.2 becomes 


f (x,y,z) 


r r f(s i* s 2 ) 


exp[2iTi[S 1 x + S 2 y + C^z^ds^ ds 2 


2.3 

2 

When S>1, (1 -C) >1 yields imaginary values for C. 
Letting = iy for real y, the term exp^iTiC^z) in equation 
2.3 becomes exp (- 2ttyz) . Waves with S > 1 are therefore 
greatly attenuated as they leave the screen and so carry no 
power away from the screen. They are commonly referred to 
as evanescent waves. For the consideration of wave fields 
at a distance of many wavelengths below the medium, the 
limits of integration in equations 2.2 and 2.3 can be 
extended to ± °° to yield 


f Cx,y,z) 




F(S 1 ,S 2 ) exp[2Tri(S 1 x + S 2 y + CjZ^ds^ ds 2 


Previously phase was measured w.r.t. the carrier at 
the diffracting medium. In evaluating the spectrum at a 
distance z from the screen, only phase measured w.r.t. the 
carrier at the point of observation will be used. This 
change is accomplished by subtracting the phase factor 2Triz 
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from the 2iriC 3 z term above to yield 

f (x,y,z) = f [ F(S 1 ,S 2 ) exp[2TTi(S 1 x+S 2 y+(C 3 -l)z)]ds 1 ds 2 

2.4 

Comparing equation 2.4 with 2.2, f(x,y,z) is composed of a 
complex angular spectrum 

G(S 1 ,S 2 ,z) = F(S 1 ,S 2 ) exp[2TTi(C 3 - l)z] 2.5 

which can be thought of as the function F(S^,S 2 ) with its 
components shifted in phase to account for propagation in 
different directions to a distance z from the screen. As 
will be shown later, it is f(x,y,z) which is actually mea- 
sured and the angular spectrum is used as a means of deter- 
mining the nature of the diffracting medium from the measured 
pattern f(x,y,z). Equation 2.4 is the expression for a 
Fourier transform pair, with z as a parameter, and can be 
represented as: 

f(x,y,z) G(S 1 ,S 2 ,z) 2.6 

If the case of a diffracting medium with no absorp- 
tion (which is a good approximation for the F-region at the 
frequencies used) is considered, f(x,y,0) will contain only 
phase variations. This is equivalent to setting A(x,y) = 1 
in equation 2.1. An attempt will be made to obtain pre- 
dictions about the properties of the diffraction pattern by 

J 
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considering a simple diffracting medium with a cosinusoidal 
variation of phase shift with position. The generalization 
of this model to more complex medium distributions will be 
then considered. The inversion of the results of this 
analysis will provide a basis for interpreting experimentally 
obtained data. 

2 . 2 Fourier Analysis of a Weakly Refracting Screen 

In the interest of simplicity only variations in the 
x direction will be considered. Let f(x,0) consist of a 
cosinusoidal phase variation imposed on a unit amplitude 
carrier given by 


f(x,0) = exp(i<£ 0 cos 2TTJ) 2.7 

where d is the periodicity of the modulation and therefore 
of the screen. If <f> Q << 1, then f(x,o) = 1 + i«J> 0 cos 2 ttj 
which can be shown to have the following angular spectrum 
(the values of S in the parentheses indicating the angle of 
the component). 

F(S) = 1 (S = 0) + — ^2. (S = + ~ 2 ~ (S = -jO 

The first term represents the carrier and the other two 
terms are phase quadrature terms. Letting the integral in 
equation 2.4 (without the y dependence) become a sum, 
f(x,z) becomes: 
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- V - 

f(x,z) = {[F(Sj) exp 2iri CC^. - l)z] exp ZttiS^x} 

j 

substituting F(S) the angular spectrum G(S^,z) becomes: 

G(S) = 1 + i<p Q exp 2iri(C - l)z 

= 1 - <j) ■ [sin 2n(C - 1) z + i cos 2u(C - l)z] 

2.8 

The first term represents the carrier. The second 
term is real and so is in phase with the carrier. It will 
therefore produce amplitude variations across the wave 
front. The imaginary term is in phase quadrature with the 
carrier and so represents phase modulations on it. The 
terms "amplitude components", and "phase components" of the 
angular spectra will be used rather loosely to refer to 
terms in phase or phase quadrature with the carrier respec- 
tively. Using this interpretation, the real and imaginary 
parts of equation 2.8 indicate that for a given value of z, 
the power transfers between phase and amplitude components 
as a function of S through the term cos 2 ttCz. This repre- 
sents the varying phase advance as a function of S in 
traveling a distance z and for d >> 1 can be expressed as 

? 1 7 

$ = 2 7T (C - 1) z = 2 ttz [ (1 - S^) 2 - 1] = -ttS z 

The phase modulation in equation 2.8 disappears for values 
of n such that irS^z = (2n - l)j, or S = [ (2n - l)/(2z)] 2 . 
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What this means is that for S << l/(2z), only phase effects 
will be present, but amplitude components will dominate the 
angular spectrum as S approaches the case for n = 1. There- 
fore for screen periodicities d near 

d = | n = 1,2, ... 2.9 

The phase modulations produced at the diffracting medium 
will appear as amplitude variations at a distance z, which 
will tend to produce nulls in any phase spectrum for the 
conditions of equation 2.9. The nulls will not be very 
sharp due to the motion of the satellite causing variations 
in the effective height of the diffracting layer with 
latitude. There is also the change in z^ , z 2 , and h Q with 
zenith angle. Both these effects tend to change the values 
of d for which the nulls occur and therefore produce shallow 
nulls. Only the first null will have much chance of being 
observable due to its large separation from the second null. 
The -case for S << 1 / ( 2 z ) will be of primary concern in any 
study of phase fluctuations because for a thin screen under 
this restriction the pattern on the ground is a simple pro- 
jection of the screen. 

Until now it has been assumed that the source is at 
an infinite distance from the diffracting screen. The 
actual case is a source at a finite distance. As was shown 
by Briggs and Parkin (1963) the change is accomplished by 
replacing z by z^z^l (z^ + z 2 ), where z^ = h^, z 2 = h Q - h^ 
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shown in figure 1. This yields the following expression 
which replaces equation 2.8: 


G(S,z) = 1 - ■ <b sin[27r(C - l)z, r-~] 

o X 


2.10 


+ id cos [ 2tt (C - l)z 1 jj— ] 

0 


For z 2 » h Q >> z^, equation 2.10 approaches 2.8 since 
z 2 / (Zi + z 2 ) -*• 1, but as z x increases w.r.t. z 2 the ratio 
z 2 /h o decreases until for z^ = z 2> z i^q ~ 1/^. This ex- 
treme case won't occur for the restricted zenith angles 

7 

encountered in our experiment but ratios on the order of ^-q 


could occur. The effect this has on equation 2.9 is to 
lower the sizes of d for which nulls occur by the factor 
Yz 2 /b 0 to produce 


d 


2 z- 


[ 7H" 




2.11 


2 . 3 Strongly Refracting Medium 

Previously only weak scattering was examined. It 
will be shown that there is also the need for considering 
cases other than 4> << 1. The consequence of this 

inequality not holding is seen most clearly by looking at 
equation 2.1. For amplitude effects (<j>(x,y) = 0) there is 
a linear relationship between f(x,y) and A(x,y) independent 
of the size of the amplitude fluctuations imposed by the 
medium. This will allow A to be broken into its Fourier 
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components and £ to be constructed by the superposition of 
the wave fields of the individual components. For a phase 
screen a linear relationship doesn't in general hold between 
f and 4> and therefore will not allow the superposition 
principle to be applied. Letting A = 1 and treating the 
one dimensional case, equation 2.1, with <j> (x) broken into 
its Fourier cosine series, becomes 

00 

f (x , 0) = exp{ i 4> n (cos } 

n=0 

00 

= exp { i [1 + ^ ^n^cospp—)] > 

n=l 

Linearity can only be assumed if E<f> n for all the screen 
components is << 1, for then the exponential can be expanded 
into the form 

00 

f(x,0) = 1 + ^ i<J> n cos 
n=l 

There is therefore a need to treat the case for <f) Q ~ 1 dif- 
ferently. As will be shown later, this case will occur 
more often for large structure sizes in the diffracting 
medium. For the larger structures the distance z 1 is such 
that simple ray optics can be applied. To a first order 
approximation, the phase pattern on the ground will be a 
direct projection of that leaving the medium. This is the 
same as for the smaller structures which obey the limit 
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<f> o << 1, and thus allows a consistent interpretation for 
both cases. 

For the sake of completeness, the case for <p Q ~ 1 
will be looked at in the same manner as that for <p Q << 1 to 
show the consequences of this limit not holding. Using 
equation 2.7 and applying standard phase modulation theory, 
f(x,0) becomes: 

00 

f(x,0) = ^ (i) n J n (<{> 0 ) exp (2Trin|) 2.12 

n = -oo 

with F(S) comprising components at S = 0, ± j, . . ., 

± j-, . . ., and amplitudes (i) n J n 0 ) • For n > the J n 

(<j> Q ) terms become very small and can be neglected, leaving 
2n significant spectral components. The effect of these 
higher order terms will be an increased weighting of the 
higher components of the power spectrum over what would be 
expected from weak diffraction theory. Due to the small 
angles of diffraction for these large structures, only the 
first few components from a particular point on the screen 
will be observed, and the components refracted from other 
parts of the screen won't reach the observer. For non- 
linear effects of this nature, the presence of more than one 
periodic pattern in electron density will not produce simply 
additive phase spectra, but will contain also the additional 
cross terms for S = m/d^ ± n/d 2 , where m, n are the orders 
of the terms in equation 2.12. 
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2 • 4 Characteristics of the Inhomogeneous F-region 

Some important characteristics of the inhomogeneous 
F-region must be considered if the phase power spectrum it 
produces are to be analyzed properly. 

Previous experimental work has determined some 
physical properties of the F-region inhomogeneities. By 
using spaced receivers, (allowing the construction of a two 
dimensional correlation function in the plane of the receiv- 
ers) the shape of the correlation contours of the irregular- 
ities were found to be ellipsoidal, with the major axis 
aligned with the local geomagnetic field. Detailed obser- 
vations of amplitude scintillation performed at Boulder, 
Colorado (49 degrees N. latitude) produced the following 
characteristics (Jespersen, Kamas 1963): 

r = transverse dimension: 0.50 km - r - 2.90 km; mean 
o o 

r = 1 km 
o 

AN q = total electron variation which increases sharply for 
northern latitudes, — 40°N. 

H = height of the irregularities; mean H = 325 km 

T = thickness of the region; T = 120 km 
This description is based on amplitude scintillation data, 
however, and biased significantly by the diffraction process 
as discussed above. 

It has been shown (Dyson 1969, Kelleher 1970) that 
the irregularities extend to heights much greater than the 
average height. It is only because the fairly concentrated 
band of irregularities near 300 km produces the greatest 
variations in electron density (peak density of the F-region) 
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0 (STATE COLLEGE, PA) 

Figure 1 Flat earth geometry of the satellite, ionosphere, 
receiver system 
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that the thin diffracting screen approximation can be used. 

It is also accepted that the axial ratio ranges from 1 to 
10, with an average of about 3. 

Some geometrical parameters will be brought in now 
for later use in Chapter 4. Because . the Beacon B satellite 
is in a polar orbit, the ray path of its signals will sweep 
along the major axis of the field aligned irregularities for 
overhead passes. At the point of observation, 40,8 degrees 
N. latitude, the magnetic field has a dip angle. of 72 de- 
grees. For overhead passes the propagation vector k will 
therefore be parallel to (geomagnetic field) , for the 
case where the satellite is at its southern most point of 
observation (zenith angle of - 23 degrees). As the satellite 
progresses northward this angle increases to a maximum of 
about 40 degrees. As found by Jespersen and Kamas (1963), 
the mean ionospheric height, thickness of the diffracting 
medium, and frequency of occurrence of irregularities all 
increase with latitude. 

2 . 5 Diffraction by an Inhomogeneous Ionosphere 

Due to the lack of amplitude information in the 
recorded data (see Chapter III), it will not be possible to 
produce an exact description of the phase power spectra 
leaving the diffraction medium, rather a qualitative 
description of the general shape of the lower end of the 
expected spectra will be sought. 

The power spectrum of the phase modulation 


17 


P (k ,k ) imposed upon the signal near a thin diffracting 
y x y 

screen and the power spectrum of electron density fluc- 
tuations p jsjC^ x >^y»^ 2 = 00 are related by 


P , (k ,k ) = 2 it ( r X) ^LP m (k ,k ,k = 0) 
d x’ y e N x’ y’ z 


2.11 


(Lovelace et al. 1970, Cronyn 1970). 

L = thickness of the diffracting region 
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r = classical electron radius = 2.83 x 10 m 
e 

X = wavelength of the signal 
kj = spatial wave number in the j direction 
k z = 0 was chosen to represent a thin diffracting 
medium. 

The shape of the two spectra are similar and thus represent 
a simple projection of the electron density variations of 
the medium on the wavefront in the form of phase flue- • 
tuations. It has been observed (Porcello and Hughes, 1968) 
that the form of P^ for the case of auroral zone inhomogene- 
ities can be given by a power law variation of the form: 


7 7 + ^ 

P M (k ,k ,k = 0) = A(k c + k *) - — 

N v x’ y* z J K x y J 2 


with 8 U found experimentally to range from 3 to 4 . 

A check on the validity of the assumption <{> 0 << 1 
will now be made. The phase deviation in traversing the 
medium is given by: 
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ra 

<j> = - r X j ANdl 
o 

where r is the classical radius of the electron. Letting 
e 

the deviation in electron density AN be Gaussian, <j>(x) be- 
comes (as shown in Briggs and Parkins 1962) 



with typical values to be expected 

ip = angle between and k = 20 degrees 

AN = maximum electron density change = 1.0% of lO 1 ^ = 

9 1 m 

10 l 3 
m 

X = 15m 

a = axial ration of the ellipsoid = 5 

r = 5 km 
o 


4»(x) = 8.3 x 10" 2 exp[-2(-^-) 2 ] 

io - 5 

This gives a maximum phase shift of <j>(0) - 0.5 radians. It 

is therefore conceivable that for larger structures, r Q = 

50 km, and larger electron density variations AN q = 10%N q , 

values of <j> could be on the order of 50 radians. Thus we 

could expect higher order components along with the first 

order term. The number of side components n is roughly 

9 1 

n = <j> . If we use the value of 10 — , for AN , the largest 

_ o 

m 

irregularity for which 4> 0 < .25 radians can be found to be 
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.25 


8.3 x 10' 2 — 2- 
10 m 


r = 3 . 0 km 
o 

If it is assumed that AN q is proportional to r. Q , there will 
be a cutoff at which the approximation 4 > q << 1 no longer 
occurs. It is therefore expected that the power contri- 
bution to the high end of the observed phase spectrum will 
contain the sum of contributions from the short periodicity 
screens and also the higher order components of those large 
structures which lie below the cutoff. 

Using equation 2.11 the amplitude effects aren't 
significant for d >> [2z^ (z 2 /h Q ) ] . For Xz^ = 300 km, 

\z - 700 km this condition becomes d >> 3 km to yield a 

id 

lower bound on structure size, above which the phase pattern 
produced at the ground is a simple projection of the pattern 
at the diffracting medium. This cutoff will be chosen to be 
5 km. 

Due to satellite motion the 2 -dimensional spatial 

diffraction pattern at the ground has a velocity V = - h^/ 

Ch - h.). Because of this motion an observer at a single 
o 1 

receiving station will observe a one dimensional temporal 
scan of this pattern, assuming that it moves without change. 
This assumption is equivalent to assuming that the dif- 
fracting region is very thin in vertical extent. The 
spatial wave numbers of the pattern k^ can be equated to the 
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temporal frequency components f ^ by k^ = 2Trfj/V or 

kj = 2 Trf j (h Q - h i )/(h i V s ) to produce k^ = 2.09 f sec/km. 

The power spectrum in frequency of the received temporal 
pattern can therefore be used to reproduce a one dimensional 
spatial pattern. 



CHAPTER III 


SIGNAL HANDLING AND DATA REDUCTION 
EQUIPMENT AND TECHNIQUES 


3 . 1 Source and Initial Signal Handling 

The source used was the S-66 Ionospheric Beacon-B 
Satellite with the following orbital parameters: 


Nodal Period 


104.8 minutes 


Inclination 


79.7 degrees 


Perigee 890 km 

Apogee 1070 km 

It produces continuous, harmonically, related 
signals at the frequencies 20, 40, 41, MHz (NASA 1963) . Use 
will be made of one of the circularly polarized modes of the 
20 MHz signal and the linearly polarized 40 MHz signal. 

The receiving and recording equipment was construct- 
ed about 8 years ago by previous workers, with the block 
diagram shown in figure 2. The data was collected over the 
period 1964-1968, and after processing was stored on mag- 
netic tape. The signals were also recorded on chart records, 
a sample of which is included in Appendix B. By using har- 
monics of a single tracking oscillator for both frequencies, 
the phase relationship between them was preserved for later 
phase comparison. The receiving antennas consist of a fixed, 
vertically directed, circularly polarized antenna and a 
fixed dipole with an east-west orientation. 

The direction of polarization of the incoming 
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Figure 2 Block diagram of the receiving equipment 










23 


signals will change with time to produce a Faraday rotation 
rate. This is caused by the systematic change in iono- 
spheric thickness and in the longitudinal component of the 
earth's magnetic field with zenith angle. For the single 
dipole antenna this produces severe amplitude variations 
which must be compensated for in the early stages of the 
electronics due to the critical signal level dependence of 
the phase detection circuitry. There will also occur phase 
discontinuities at the times of cross -polarization between 
the incoming wave and the antenna. These discontinuities 
will produce spikes in the temporal . diffract ion pattern 
and therefore produce a spurious addition to the phase 
spectra. Because the 20 MHz circularly polarized antenna 
was fixed it will appear to be ellipitically polarized for 
oblique incidence of the signal. This will produce an 
apparent, non-uniform Faraday rotation rate which in turn 
produces phase modulations. The magnitude of this effect is 
increased with increasing zenith angle which makes necessary 
the limiting of satellite longitudes to the range 75-80 
degrees west and of latitudes to the range 35.5-46.5 de- 
grees north. Since the above effects are reduced by a 
decreased Faraday rotation rate which is dependent on the 
integrated electron density of the ionosphere, only night 
time records will be used, 

A brief note will be added here as to why linear, 
rather than circular polar ization . was used for the 40 MHz 
signal. Both 20 and 40 MHz circularly polarized signals 



were stored on the same tape channel. There was some dis- 
tortion in the receiving equipment which generated a second 
harmonic of the 250 Hz stored signal. This harmonic was 
inseparable from the 500 Hz signal, and caused spurious 
phase modulation of it. 

3 . 2 Analog Processing 

The approach used is to compare the frequencies of 
the two previously stored signals yielding a continuous 
temporal frequency record. This represents the derivative 
of the temporal phase function w.r.t. time recorded at the 
receiving station. The frequency record represents the 
deviation w.r.t. time of the signals from their harmonic 
relationship, which will give the phase path difference 
experienced by the signals in traversing any inhomogeneities 
in the F-region. Utilization of harmonically related sig- 
nals allows the elimination of tracking errors, source 
stability errors, and phase noise induced during the record- 
ing and playback process. This is accomplished by frequency 
doubling the 250 Hz signal and comparing it with the 500 Hz 
signal to eliminate any modulation proportional to frequency, 
of which all the above produce. The 500 Hz signal was off- 
set by 180 Hz (chosen by trial and error to produce the 
least harmonic distortion) to facilitate phase detection. A 
block diagram of the electronics is shown in figure 3, where 
C.f. and B.W. are abreviations for center frequency and band- 
width respectfully. The following requirements were placed 
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Figure 3 Block diagram of filtering and discriminator circuitry 
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on the electronics to account for physical conditions and, 
or previous recording methods. 

Because the frequency discriminator used was 
sensitive to amplitude variations, severe amplitude limiting 
was imposed prior to frequency detection. This need was 
made more demanding by the deep nulls in signal caused by 
Faraday rotation, especially on the 500 Hz signal. 

Predetection bandwidth requirements were imposed in 
early limiting to prevent harmonic distortion from reaching 
the frequency detector. Minimum bandwidth (B.W.) require- 
ments are set by the phase modulation induced on the signal. 
Due to the changing effective thickness (shown in figures 1 
and 4) of the ionospheric layer with satellite latitude B, 
we will have a linear ramp imposed on the discriminator out- 
put. Hereafter this effect will be referred to as the 
apparent doppler dispersion. Its magnitude is approximated 
by: 


<J> = -r X / N(x)dx 
e o 


f = 




V d 


“ art 


(N o Jl) 


3.1 

3.2 


3.3 


N q = average electron density (assumed here to be constant) 
i = arbitrary ionosphere thickness 



[Ton (V s £/l8ho)] vs. £ 


f(£) vs. /3 


f(/9) = Sm/j[cos£(r 2 +r 2 - rr, Cos/?)-rr,j[l- 

(Sin 2 /?/ r 2 +r 2 - 2 r r, Cos/3)j 


Frequency modulation caused by the apparent doppler 

dispersion for (a) flat earth, (b) curved earth 

geometry. R = radius of the earth, r ■ (R + h.)/ 

(R + h , r. e = (R + h.)/R , * maximum 6 e of 1 

„ve o ? 1 e i e’ M 

observation 
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and becomes 


r d 

£ = • 4 AN V 4 sin 6 3.4 

2ir oh 

o 

The records were restricted to values of 6 — 6 = 

m 

23°. Using nighttime values of N Q given by Solomon (1965), 
the maximum frequency deviation for the frequency doubled 
250 Hz signal at the extreme zenith angles of the data was 
found to be f max - ± 0.5 Hz for the 20 MHz signal. Because 
both the 20 and 40 MHz signals are modulated by the iono- 
sphere similarly but with different magnitudes, A£q - A^q 

should have been used rather than A» n which reduces f to 

20 max 

about ± 0.4 Hz. 

The tracking errors were measured by previous 
workers to be as large as 10 Hz and must be considered in 
B.W. considerations of the filtering prior to the mixing of 
the two signals. Because the analysis treats ionospheric 
induced frequency modulations up to 3 Hz, (actually up to 
8 Hz, but this region wasn't analyzed because of severe 
attenuation due to postdetection filtering) this will have 
to be added into both predetection and postdetection B.W. 
limits. Combining the above, the predetection B.W. needed 
is 2 (10+3+0 . 5)Hz = + 27Hz. It was set to 30 Hz. Post- 
detection integration contained a cutoff at 2 Hz. 

The discriminator output represents the frequency 
difference of the 2 signals as a function of time. This 
necessitates that the delay of each signal through the 
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electronics be identical, i.e., that the rates of change of 
phase shift vs. frequency for each signal channel be equal. 
This was achieved with an additional filter in the 250 Hz 
channel to yield slopes differing by less than 5%, repre- 
senting a phase error of ten degrees for a frequency change 
of 28 Hz. 

Because the discriminator used was very sensitive to 
amplitude variations, severe limiting was used and signal 
amplitude was monitored to insure against fluctuations. 
Discriminator output remained constant for the following 
ranges of signal input: 

Amplitude 250 Hz: 0.6 - V(250) - 3.2 VP-P 
Amplitude 500 Hz: 0.2 - V(500) - 3.2 VP-P 
The discriminator output was placed through a 3 stage, low- 
pass, filter with cutoff at 2.0 Hz. 

Its response was linear over the range used (± 5 Hz) 
with a sensitivity of 2^-- and a noise level (due primarily 
to intermodulation in the earlier electronics) of ± .02 Hz. 

A sample of a typical chart record produced by the 
above is shown in figure 5. This information consisted of: 
Channel 1: WWV 

Channel 2: discriminator output f(t) in Hz 

Channel 3: discriminator input amplitude monitoring 

Channel 4: 40L amplitude monitoring 

3 . 3 Data Analysis 


The aim of this stage of analysis is to obtain a 
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Figure 5 Sample Record of Discriminator Output f(t) 
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power spectrum of the phase variations imposed on the 
received signal. The method used is outlined by Blackman 
and Tukey. 

To avoid the complexity required to obtain the 
frequency . components by analog methods, the data will be 
digitized.. This action replaces the continuous function 
f(t) by f(tj), where t^ represents the discrete time values 
j At , j = 1,2 . . . N, NAt = length of the record. As noted 
earlier, there exists a D.C. trend on the data. To elimi- 
nate this large, low frequency contribution, the best linear 
fit line was subtracted from the data point by point. De- 
termination of the line defined by the quantity (a + BjAt) 
was achieved by the use of : 

N 

^ [f(t.) - a - Bj] = 0 3.5 

j=l 

N 

£ CJ - f) [fCtj) - a - Bj] = 0 3.6 

j=l 

Solving for a and B, the trend and mean can be subtracted 
out. The autocorrelation function of the detrended data 
was found, and its Fourier transform taken, yielding the 
power spectrum in frequency, which will be shown later to be 
related to the power spectrum in phase. Several aspects of 
digitizing data must be considered before the usefulness 
of the data can be assured. 

Computation was performed on a PDP-10 digital com- 
puter with a limit of 23 k user available core. This set 
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a practical upper limit on the size of the data array and 
therefore the sampling rate. Because the maximum observable 
frequency with a sampling rate of N per second is 1/(2N), N 
was set so 1/(2N)=8 Hz. The data was already limited to 
~ 3 Hz due to final filtering, so no loss of information 
occurred. N therefore is equal to 16 samples per second. 
Previous work (Porcello and Hughes, 1968) has shown the 
power spectra to fall off with increasing frequency which 
will further reduce aliasing. The longest lags obtained are 
determined by the length of the record, and are usually 
taken to be 5-10% of the record length. We chose lags up to 
32 seconds which is - 16% the record length. This gain in 
range of data is seen in equation 3.7 to reduce the stabil- 
ity of the results (represented by K) . 

Maximum aliasing is determined by the attenuation of 
the output filter. The power at the folding frequency of 
8 Hz is compared with that at 0 Hz and was found to be down 
68 db . This is more than sufficient to eliminate any 
aliasing problems. 

Due to changing ionospheric conditions between 
orbits, nons tationarity between records is assumed, which 
forces statistical independence of the records, rather than 
an average over an ensemble. For the record length avail- 
able and maximum lag used, the number of degrees of freedom, 
K, is found to be (Blackman and Tukey) : 

. 2 V . 2 t T n - W3TJ 
“T“ f 


K 


3.7 



33 


T n = total record length 

T = maximum lag desired 
m 

K = 11-1/3. Use of Table II in Blackman and Tukey leads to 
an 801 confidence limit of 4.8 db. This value is only 
characteristic of the runs taken, with exact values given 
later for individual records. 

Past evidence (Porcello and Hughes, 1968) has shown 
the phase spectral density to follow a power law relation- 
ship, given by 

e H - 1 

VW “ (k x 2 + k y 2) ^ 

with k^, ky the wave numbers and B H an experimentally 
determined variable. This indicates the need for pre- 
whitening. The data being originally intended for other 
uses, prewhitening wasn't performed. This will reduce the 
signal to noise ratio. at the high end. of the spectrum, 
making that region more susceptible to interference or noise. 

The process of taking a finite sample is equivalent 
to multiplication of the infinite data set with a square 
pulse. This will accentuate the high-frequency end of the 
spectra. This effect was reduced by passing the auto- 
correlation function C(T^) through a Hamming spectral window 
D(Tj) where 

D(T.) = .54 + .46 cos<Vr./T ) | T . I < T 

j K J m' 1 j 1 - m 


0 
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The attenuation due to final integration was compensated for 
by direct measurement of the characteristics of the low-pass 
filter and multiplication of the frequency power spectra 
lines by the appropriate value. 

The flow diagram and computer program are included 
in Appendix A- 2. 



CHAPTER IV 


EXPERIMENTAL RESULTS 

4 . 1 Experimental Details 

The individual orbital conditions and statistical 
parameters of the data used are given in Table I. Due to 
the severe limitations on longitude and the restriction to 
nighttime records, only a small fraction of the available 
library of records was considered of which only 8 satellite 
passes produced significant frequency variations observed at 
the discriminator output. The effectiveness of the trend 

removal performed on the digitized discriminator output was 

6 7 

about a 10 reduction in mean value, and a 10 reduction in 

slope. Because both the original mean and trend were of the 
same order of magnitude as the desired information, these 
reductions brought their effect well below the noise level. 
For an. 80% confidence limit, there was typically a 5 db un- 
certainty range in the relative power spectra. 

A sample of the discriminator output is shown in 
figure. 5. The 40 MHz amplitude channel was displayed so 
as to observe the magnitude of the frequency spikes produced 
by the Faraday rotation and any possible amplitude effects 
it could cause. These finite width spikes have a minimal 
overall effect on the spectra due to their infrequent occur- 
rence. The scale is 1.0 Hz per division on the discrim- 
inator channel. 
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4 . 2 Spectral Density of the Frequency Modulation 

(a) Analysis of the high end of the spectra. 

Sample plots are shown on log-log scales in figures 
6-11. Figures 6 and 7 show a spectrum before smoothing and 
after smoothing. Smoothing was accomplished by taking a 
running average over three unsmoothed logarithmic spectral 
values, which themselves were arithmetic averages over two 
adjacent spectral values. The remainder of the figures are 
for smoothed plots only. This doesn't reduce the useful 
information because as is seen in figure 6 the features pre 
sent in the unsmoothed version and not in the smoothed 
version are within the range of uncertainty imposed by the 
80% confidence limit shown on each graph. The consistent 
appearance of nulls at the high end of the spectra is a 
feature to which attention is now drawn. 

For a source at a finite distance from the observer 
it was predicted in equation 2.10 that nulls in the power 
spectra of structure sizes should occur for dimensions 
given by equation 2.11. A conversion from spatial to 
frequency spectra produce nulls for frequencies 


f = 


V h. 
s 1 

“d K7 


V h. 
s 1 

OcT 


V 


r 2n ■ 

[ T 


•K“Ur 

0 0 


h. 

1 

- h.) 
i J 


) ] 2 n = 


1 , 2 , . . 

4.1 


This produces nulls at the frequencies (for n = 1,2,3, . . 
respectively) f = 0.675, 1.31, 1.69, . . . Hz. Figures 
6-11 show definite nulls in this frequency range, but 
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Figure 10 Smoothed power spectrum for orbit Figure 11 Smoothed power spectrum for orbit 
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generally not at these precise values. 

Upon a closer analysis of the geometry involved, it 
becomes clear why the predictions above aren’t realized. 

The values of the predicted nulls depend on the height of 
the diffraction medium and the distance to the satellite. 

As shown in section 2.4 the medium height increases with 
increasing latitude and h Q increases with zenith angle. The 
finite thickness of the diffracting layer also produces 
inhomogeneities at different heights. Assume for example 
that lu varies over the range 300-400 km. This produces a 
first order null ranging from 0.675 to 0.504 Hz to produce 
a broad null centered about 0.59 Hz. It could also be the 
case that there exist diffracting layers at different 
heights. This effect could come about through a patch 
structure of inhomogeneities as observed by Kelleher and 
Sinclair (1959). Here structures appear in groups pro- 
ducing the possibility of patches occuring at different 
altitudes, 300 and 400 km for example. This would also 
produce distinct first order nulls at 0.675 and 0.504 Hz. 

No quantitative analysis will be made of this region 
for the following reasons: (1) the confidence limits are of 

the same order of magnitude as the depth of the nulls, 

(2) the data is believed to be approaching the system 
noise level in this region as indicated by the character- 
istic flattening of the spectra. Had these two conditions . 
not been present, some estimate of the height and thickness 
of the diffracting medium could possibly have been obtained 
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from the position and. shape. of the nulls. It is interesting 
to note that although. the nulls should occur in both 
amplitude and phase power spectra, amplitude studies 
(Rufenach 1972) have only rarely observed. them, while they 
are quite evident in this study using. phase scintillation. 

(b) Analysis of the low end of the spectra. 

This is the region of the spectra containing fre- 
quencies below the first null. This frequency is found, 
using equation 4.1 for n = 1 and converting from normalized 
distances to actual distances, to be 


f = 


£ 


h 


X] = V, 




"1 - 

rnr- 

2 O 


4 . la 


It is in this region that some insight might be gained as to 
the physical nature of the diffracting medium. 

From the records studied there tends to be a definite 
slope 3 in the log-log plots shown. This would indicate a 
frequency power spectrum of the form 


P(f) - Cf)' B 


4.2 


Upon taking a visual estimate of the slope, it was found 
to have a mean value of 0.71 and standard deviation of 0.10. 
These values were found neglecting the 3 value for orbit 
10566 due to its large deviation from the norm. Upon look- 
ing at the chart record for this run a considerably larger 
amount of low frequency variations were present, thus 
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indicating a greater than average (taken over the 8 passes 
used) number of larger structure sizes w.r.t. the smaller 
sizes . 

(c) Statistical stationarity of the results. 

Due to the similarity of values of 8 for the dif- 
ferent orbits, some support is rendered to the case for 
statistical stationarity of the medium w.r.t. this parameter. 
Porcello and Hughes (1968) also. found similar values for 8 
for structure sizes smaller than considered here. From the 
limited sample taken, there is the indication that the 
parameter is a characteristic of the medium and not depen- 
dent upon temporal changes occur ing in the inhomogeneous 
makeup of the diffracting medium. Orbit No. 10566 produced 
a 8 value with considerable deviation from the mean value. 
Because of the limited sample taken, no conclusions will be 
drawn from this deviation. 

4 „ 3 Discussion of the results 

It was shown from the graphs and expressed in 
equation 4.2 that the power spectrum obeys a power law 
variation with frequency. What is of more physical signi- 
ficance is the power spectrum of the phase difference c{> ( t ) 
of the two signals. Noting that f(t) = d<j>(t)/dt, the 
Fourier transform of f(t) is 

00 

Ff(f) = J f (t)e' 2TTl£t dt 
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F f( £ ) = f 3 * t 4>(t)e" 27rl£t dt 

- 00 

= i f/“<K.t)e- 2,,ift dt 
- 00 

= - ifF^Cf) 4.3 

the comparison of power spectrum yields: 

P f tf) - |F f (f)| 2 - f 2 |F^(f)| 2 = f 2 P + (f) 

P^Cf) - f' 2 P f (f) 4.4 

P^tf) « f - C6* 2 > 4.5 

It is now desirable to determine the nature of the spatial 
power spectrum in phase leaving the diffracting screen, 
denoted by P^'(k). Equation 2.8 shows how the angular 
spectrum in phase changes as the distance from the screen 
increases. For the cosine diffracting screen considered 
previously, the power at the screen is given by 

P^CS) - | F (S) | 2 - 1 ♦ * 0 2 (S = £) 

On the ground the power contained in the phase part of the 
angular spectrum is 

P^'CS) = | G (S) | 2 = 4> o 2 cos 2 [2Tr(l - 


S 2 1 

*7 - i) 2 ] Cs = J) 
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P. ' (S) = <f> 2 cos 2 uS 2 z 4.6 

<p o 

Expanding the cosine term P^'(S) becomes 

|P^'(S)| - 4 0 2 (1 - ttzS 2 ) 4.7 

2 

Because it has already been assumed that zS tt << 1, the 
spectrum is independent of S and identical to that for 
z = 0. Next the dependent variable in equation 4.5 will be 
changed from frequency to wave number along the direction of 
motion k using k « f 

P^ (k) « k' (B+2) 4.8 

A proportionality rather than an equality is used due to the 
nature of measurement of the power spectra. The power 
spectrum of the electron density fluctuations is therefore 

P N ' (k) « k" (e+2) 4.9 

What is obtained therefore is a power law variation 
in electron density fluctuation over the range of iono- 
spheric structure sizes 5.0 km to 73 km. The character- 
istic exponent was - 2.71. Previous work by Rufenach (1972) 
found using amplitude scintillations, that 8 ranged from 
0 to 1 for structure sizes 0.6 to 4 km. Using phase 
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variations Porcello and Hughes (1968) found 8 for 2 records 
to be 0.8 and 1.0 for auroral conditions. 



CHAPTER V 


SUMMARY 


5 . 1 Conclusion 

A study was made of the inhomogeneous nature of the 
nocturnal F-region from a ground station at 40.8° N. lati- 
tude using phase information from two harmonically related 
satellite signals. A phase comparison of the two signals 
eliminated tracking and processing errors common to both, 
but not the dispersive ionospheric phase modulation. From 
the phase modulation of the signal, the power spectra of 
the electron density variations in the diffracting medium 
(of wavenumber k) was found to follow a power law relation- 
ship P N '(k) = k‘ ( - 6+2) . Values of 6 were found to lie in the 
range 0.50 to 1.5 with a mean of 0.71 and standard deviation 
of 0.10. Because of experimental constraints, this expres- 
sion could only be verified over the range of structure 
sizes ( 2 it / k ) from 5.0 to 7 3 Km. 

What is believed to be Fresnel oscillations was also 
observed in the diffraction pattern. The characteristic 
nulls in the power spectra of phase variations were present 
at near the expected frequencies, but due to experimental 
uncertainties no detailed analysis was performed on them. 

Statistical stationarity of the medium w.r.t. the 
parameter 8 was considered., Because of the small standard 
deviation of this value, the exponent - (8 + 2) appears from 
the limited sample of ionospheric conditions observed to be 
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a characteristic of the inhomogeneous F-region. One orbital 
pass did yield a value for 8 of 1.5, but further measure- 
ments would be necessary to examine its significance rela- 
tive to the above stationarity . 

The advantage of the use of phase over amplitude 
information was made clear by the direct interpretation 
given here to the data in terms of parameters of the dif- 
fracting medium. Amplitude studies have the problem of 
interpreting the effects of diffraction, as do phase studies 
for structure sizes below a cutoff value. It is the larger 
structures which allow observation in the near field and 
are therefore without major diffraction effects. 

5.2 Suggestions for Further Work 

This study of the spatial power spectrum of electron 
density fluctuations in the ionosphere has been constrained 
in several ways by experimental limitations, and the method 
could be improved by easing these conditions. 

For the smaller scales, less than a few kilometers, 
diffraction effects complicated the interpretation of data. 
The use of higher frequencies would extend this range and 
would also reduce the depth of phase modulation in the iono- 
sphere and thereby lead to a more direct relationship be- 
tween the diffraction pattern and the medium. However, 
unless the system noise level could also be reduced for 
example by the use of higher powers or higher gain (tracking) 
antennas, and improved recording methods, the extended 
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range might well be lost in noise due to the lower sen- 
sitivity of the effect on higher wave frequencies. 

For the large scale end of the spectrum the study 
has been limited by the finite length of record which can 
be obtained in the limited field of view of a receiving 
station. Although some small extension of the range 
(perhaps a factor of two) might be possible through the use 
of tracking antennas, this would introduce greater changes 
in the experimental geometry and further complicate the 
analysis . 
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APPENDIX A 


1. The following program was used to take a con- 
tinuous signal and sample at a uniform rate to produce a 
data array for digital analysis. The sampling was control- 
led by a D.C. level on a separate channel of the tape re- 
corder. The comparator (number 34 referenced in lines 
170,200 of the program) sensed when the level was + 1 or 0 
VDC to control the interrupt subroutines. The comparator 
was fed into an ’AND' gate along with the clock set at 16 
pulses per second. The output was then fed to interrupt 

0 (referred to in line 80) . The data was fed into the A/D 
converter channel 0 (referred to in line 230). The data was 
stored on Dectapes under the file names TEMPI . DAT , TEMP 2 . DAT 
etc. in arrays A(4100). 

2. The flow diagram in figures 12, 33 is intended 
as an outline of the method used in data handling and not 
for use in detailed programming. The following Scientific 
Subroutine Package programs were used. 

Auto: calculates the autocovarience of the data array A 

RHARM: calculates the Fourier coefficients of the real 

data array A. 

6X 
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REAL A( 4100 ) 

LOGICAL RSCS 
INTEGER A2(8) 

INTEGER A1C4100) 

EQUIVALENCE (A,A1) 

EXTERNAL ICRAC 
DO 100 1=1,8 
CALL HI NIT 

CALL GPI (0, I ERR, I CRAC , 0, 4 100 , A1 ) 

CALL STINT 

C START AN ALOG: SET MODE CONTROL TO OP, RUN 

CALL SAMO('OP') 

TYPE 1000 

1000 FORMAT ( 'START RECORDER', $) 

ACCEPT 1001, IANS 

1001 FORMAT ( I ) 

CALL SLMOC ' RUN ' ) 

105 I F( RSCS ( 034 ) ) GO TO 110 

GO TO 105 
110 CONTINUE 

I F ( RSCS ( 034 ) ) GO TO 110 

CALL GP I ST A (0, .TRUE,, RETURN, 41 00, INDEX) 

TYPE 500 

500 FORMAT (' DECTA PE NO. FOR OUTPUT?', $) 

ACCEPT 501, NOT 

501 FORMAT ( I ) 

NDT=NDT +8 

GO TO (1,2, 3, 4, 5, 6, 7, 8), I 

1 CALL OFILECNDT, 'TEMPI') 

GO TO 9 

2 CALL OF I LE ( NOT, ' TEMP2 ' ) 

GO TO 9 

3 CALL OFILECNDT, 'TEMP3' ) 

GO TO 9 

4 CALL OFILECNDT, 'TEMP4') 

GO TO 9 

5 CALL OF I LECNDT, 'TEMP5 ' ) 

GO TO 9 

6 CALL OFILECNDT, 'TEMPO ' ) 

GO TO 9 

7 CALL OFILECNDT, 'TEMP7') 

GO TO 9 

8 CALL OFILECNDT, 'TEMP8') 

9 WRITECNDT, 10) INDEX, ( A 1 ( J ) , J=l, INDEX) 

10 FORMAT C ' ' , T40 , 15/' ®,20(1X, 15)) 

A2CI ) = I NDEX 

100 CONTINUE 



54 - 


REAL A(4100) ,R(512) ,COEFIC(50) 
DATA (COEFIC ( J) , J=1 , 50/50 VALUES/ 



FIGURE 12 


FLOW DIAGRAM OF DATA ANALYSIS PROGRAM 
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APPENDIX B 



Figure 14 Sample of a record of the processed data which 
is stored before processing on magnetic tape 
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I ) CONTINUED 


I — — 

L=NUMBER OF LAGS DESIRED=512 
TM=MAXIMUM LAG IN SEC. =32.0 
T=SAMPLE PERIOD IN SEC.=?0. 0625 
M=IMPUT PARAMETER FOR RHARM=8 



FIGURE 13 FLOW DIAGRAM (CONTINUED) 










